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Abstract 

The influence of surface tension and size on rounding of cell aggregates 
are studied using chick embryonic cells and numerical simulations based 
on the cellular Potts model. Our results show exponential relaxation in 
both cases as verified in previous studies using 2D Hydra cell aggregates. 
The relaxation time decreases with higher surface tension as expected 
from hydrodynamics laws. However, it increases faster than linearly with 
aggregate size. The results provide an additional support to the validity 
of the cellular Potts model for non-equilibrium situations and indicate 
that aggregate shape relaxation is not governed by the hydrodynamics of 
viscous liquids. 

1 Introduction 

It is now well established that certain living tissues behave as viscoelastic fluids 
PP. Classical experiments on the relaxation of embryonic aggregates subjected 
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to mechanical deformations, either by centrifugation |2] or by compression [3J, 
demonstrated that tissues relax as elastic materials on short time scales and as 
viscous liquids on long time scales. The relaxation of compressed cellular aggre- 
gates was analyzed using a Kelvin model of viscoelasticity [3J . As stated by the 
authors, this model may represent a strong oversimplification of the tissue rhe- 
ology but it gives an estimation of the tissue elastic and viscous parameters |3J. 
The fusion or the rounding up of cell aggregates are similar to the coalescence 
or rounding up of liquid droplets. Gordon et al. showed that these processes 
are driven by tissue surface tensions a and opposed by tissue viscosity 77 mak- 
ing a formal analogy with the rheology of viscous liquids. They obtained first 
estimates of the ratio rj/a for chicken embryonic cells However, the authors 
gave only orders of magnitude of this ratio as they performed only preliminary 
experiments, presumably because it was difficult to obtain suspended three di- 
mensional aggregates which were fixed over long times (rotation or flattening on 
a surface are undesirable). Later, the surface tensions of different chicken em- 
bryonic tissues have been measured using a parallel plate compression apparatus 

In previous investigations one of us showed that a two dimensional (2D) 
assembly of Hydra cells behaves like a fluid: (i) cells exhibit collective [7| |S] or 
random [5j motion, (ii) the hydrodynamical laws also apply for the rounding-up 
of initially elliptical aggregates |S] • The shape relaxation was analyzed using the 
2D hydrodynamical equations derived by Mann et al. for polymer monolayers in 
the case of dissipation dominated by surface viscosity ^UJ. This method allowed 
an estimation rj/a ~ 0.5 min/^m although a was not measured independently. 
Compared to the relaxation analysis of free 3D aggregates, 2D aggregates offer 
some advantages: they are not moving, nor rotating or flattening on a nearby 
surface. The 2D system also allows to track simultaneously cell motion and a 
direct comparison with 2D simulations. 

In this paper, we summarize the main results obtained previously on 2D 
Hydra aggregates [Sj. We also present current experiments on the rounding 
of 3D aggregates of chicken embryonic cells (neural cells from the retina) and 
simulations of rounding of 2D aggregates using the cellular Potts model. Because 
of its flexibility and simplicity of implementation, Potts model has become a 
common technique for cell level simulation of biological tissues. It describes 
effectively the global features of tissue rearrangement experiments, including cell 
sorting and even the cell movements during the entire life cycle of the slime 
mould Dictyostelium discoideum |12| . In this study, we show that Potts model 
effectively captures the exponential relaxation dynamics of aggregate rounding. 
However, both new experiments and simulations show a deviation from the 
hydrodynamics of viscous liquids. 
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2 Results 



2.1 Rounding of 2D Hydra aggregates 

Hydra aggregates were prepared according to previously reported methods 7 
with 25 \xm spacers (Fig. 1A). A sequence of events in a pure endodermal ag- 
gregate initially elliptical is displayed in Figs. 1B-C. Rounding up is completed 
after 6h. We studied also the rounding of mixed and pure ectodermal aggregates 
Images were digitized at intervals of 5 ~ 10 min and analyzed on a computer 
using the software NIH Image ^3] ■ The deformation parameter d = — — 1 was 
used to describe the aggregate shape, where M and m are the major and minor 
axis of the fitted ellipse, d is zero for a perfect circle and increases for elongated 
aggregates. 

Fig. ID shows a log-linear plot of the time dependence of the d for the differ- 
ent 2D aggregates investigated. For each aggregate, d decreases exponentially 
with time. Such exponential behavior is expected theoretically in 2D from the 
hydrodynamics laws in case of viscosity dominated by the surface viscosity and 
rounding driven by line tension . In case of 2D aggregates confined between 
two glass slides, if the cell/glass friction may be neglected compared to cell/cell 
friction, the relaxation time is given by relation: 

T c = %R , (1) 
a 

where R is the radius of the rounded aggregate (final time), 77 is a bulk surface 
viscosity, and a the tissue/external medium surface tension [8J. The typical time 
scale of rounding depends thus on the ratio between the surface tension and 
viscosity. The aggregates round to decrease superficial tension, while viscosity 
acts in the sense of delaying the rounding. Fig. IE shows that for each aggregate 
type, T c is increasing with R as expected from the hydrodynamical law Eq. ^ 
The ratio 77/(7 = 0.5 ± 0.2 min//im which gives r\ = 3 x 10 4 Pa.s (1 Pa.s = 10 
Poise) using <r=l mN/m as a typical value of tissue surface tension [3J[51I§]- 
This is 10 7 times higher than the water viscosity. 

2.2 Rounding of 3D chicken embryonic neural retina ag- 
gregates 

We performed relaxation experiments on chicken embryonic tissues which have 
been already characterized by surface tensiometry [5]. Neural cell tissues 
were obtained from the retina of 9 day old chicken embryos. After dissection 
and trypsinization, cells were reaggregated by shaking them slowly in 12-well 
plates (200 r.p.m.) containing culture medium (BME culture medium with 
10% fetal Calf Serum [TJ) at 37°C during 16 hours. We could not obtain 
proper 2D chicken aggregates for relaxation experiments. We used instead non- 
compressed 3D aggregates, just floating near a bottom slide glass (Fig. 2A). 
These aggregates are sufficiently stable for shape analysis, they do not move or 
rotate significantly and do not spread on the bottom slide glass. Figures 2B-C 
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show the shape of a 9-day old 3D aggregate between initial and final time (24h) . 
The time dependence of the deformation parameter of these aggregates is well 
exponential as expected from hydrodynamics with a relaxation time T c of about 
1 day (Fig. 2D). In case of rounding of 3D ellipsoid droplets, T c is given by a 
slightly different relation than the 2D equation ^ 4 : 

T c = (3^R , (2) 
a 

where the constant j3 ~ 0.95 in case of aggregate internal viscosity much larger 
than external medium viscosity [3], and R is the radius of the rounded 3D 
aggregate. 

At a first sight, T c seems to increase linearly with aggregate radius. The 
doted line in Fig. 2E corresponds to such a linear fit with Eq. (correlation 
coefficient, 0.880). The ratio rj/a is 16.8 ± 4.5 min//im. The surface tension 
of neural aggregates dissected at 6 days and maintained in culture for various 
time periods was measured with the compression plate apparatus [H] . The value 
increases with continued incubation for about a day (between 6 and 7 days), 
then stabilizes for a time period of about 2 days at a — 1.6 ± 0.6 mN/m and 
increases again between 8 and 9 days to about a = 4.0 ± 1.0 mN/m. Taking 
this latter value, we obtain the following neural tissue viscosity value at 9 days: 
r\ ~ (4 ± 2) x 10 6 Pa.s. It is much higher than in Hydra aggregates. This is 
not surprising. First, because the kinetics of rounding is much slower in chicken 
embryonic neural aggregates. Also because cell packing is much higher in these 
3D aggregates (Fig. 2B-C) than in the 2D Hydra aggregates (Fig. 1B-C). It is 
likely that tissue surface tension and tissue viscosity depend both on cell/cell 
adhesion [5j and that adhesion is higher for 3D chicken embryonic aggregates 
than 2D Hydra aggregate. 

However, T c is probably increasing faster than linearly, i.e., a power law 
~ R 152 fits better the data (correlation coefficient, 0.939). Although there is 
large dispersion of the experimental points which is often inherent to such a 
biological system, the visual impression is also better than with linear fit (Fig. 
2E) . Of course, these data should be confirmed but they seem to indicate that 3D 
chick embryonic cell aggregates do not behave as viscous liquids. The relaxation 
of initially elongated aggregates toward a circular shape is clearly exponential 
but relaxation time is not proportional to the aggregate size. 

2.3 Simulations 

The model is based on the cellular Potts model and is defined as follows: at 
each site of a square lattice with dimensions L x L we attribute a spin that 
may assume any integer value. The set of all sites with equal spin S defines a 
cell, labeled by S. 

Adhesion between cells originates in the interaction between each spin and its 
20 surrounding neighbors (up to fourth neighbors to avoid pinning to the square 
lattice The lowest interaction energy, here taken as zero, happens between 

equal spins, simulating the absence of surface tension between sites belonging 
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to the same cell. Energy between different cells is taken to be positive and 
when different types of cells are considered, the interaction energy may also be 
different, depending on the types involved in the interaction. The intensity of 
each possible interaction must be specified. Hence, for example, for a system 
made of an aggregate of cells of two types only and medium we must define three 
parameters. The result is a proportional to the shared interface length, with 
a proportionality constant depending on the neighboring cell types to reflect 
differential adhesion strengths. 

Cell resistance to compression is modeled by an energy term proportional to 
the square deviation of cell area from its target area At{S). The proportionality 
parameter A plays the role of a Lagrange multiplier for this constraint and 
regulates cell compressibility. The complete energy used in the Monte Carlo 
protocol is 

H = EE^5 iV (l - S 8ii ,s tli> ) + A]T [a(S) A T (S)} 2 , (3) 

ij i'j' S 

where a(S) is the area of cell S 1 , Sij is the spin at site and Eg^^^-, is 

the interaction energy between neighboring sites labeled by Sij and Si'ji . If 
Syj = Si'ji then Eg^ s.,., = since in this case both sites belong to the same 
cell. 

The aggregate surface tension (cr, at T=0) derives from the coupling energies 
in Eq. by [TT] 

a = E cM -^, (4) 

where E c m and E cc are the adhesion energies per unit of interface between cells 
and medium and between cells, respectively. 

We implement the simulations as follows. We randomly choose a lattice site 
and one of its eight surrounding neighbors. We then propose that the site 
assumes the neighboring spin value S' and accept with probability P: 

( exp(-Ai?/T) if AH > , 
P(S ->S') = < 

{ 1 if AH < , 

where AH is the variation in energy produced by the change of spin values. 
A Monte Carlo Step (MCS), is defined as L 2 exchange attempts, the total 
number of sites in the lattice. For sake of simplicity, we use the MCS as the 
unit of "time" in what follows. In fact, the random sampling of spins, and the 
immediate dissipation of their energy (perfectly local, both in time and space) 
imply that the correspondence with real time is not straightforward 

In this model the simulation temperature T relates to active cell membrane 
motion. Higher T implies higher membrane activity and vice-versa. T simulates 
membrane fluctuations driven by the cytoskeleton in real cells and should not 
be confused with much smaller thermal fluctuations 

Aggregate rounding was investigated using as initial states elliptical aggre- 
gates with different numbers of cells: 544, 953, 1642, 2204, and 3233 cells. Cells 
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areas have approximately 100 lattice sites and cell-medium surface tensions have 
values: 4.5, 5.5, 6.5, 7.5, 8.5, 9.5. To generate these surface tensions we keep con- 
stant the interaction energy between cells, E cc —1, and vary the energy between 
cells and medium, E c m = 5, 6, 7, 8, 9, and 10. The simulation temperature 
used is 10 and A = 1. For aggregates with 1642 cells we simulated all of the 
above surface tensions, however for other aggregate sizes we studied only the 
three lower values of surface tension. In Figs. 3A and 3B we show, respectively, 
the initial and final state (complete rounding-up) of the simulated aggregate 
with 544 cells. 

We calculate the deformation in the simulated aggregates using the ratio 
of the eigenvalues of the two dimensional variance of mass distribution of the 
aggregate that is defined as 



where x and y are the coordinates of all lattice sites belonging to cells in the 
aggregate. As the aggregate rounds, the ratio between the two eigenvalues of 
the matrix above, e\ and e-i goes to 1. In analogy with experimental studies, 
we define aggregate deformation using the ratio between the larger and lower 
eigenvalues of the matrix minus one: d — e\jei — 1, (ei > e-i). This measure of 
the deformation differs from that used in the experiments but should yield the 
same exponential decay. 

In Fig. 4 we plot d as function of time for aggregates with 1642 cells for 
surface tensions 4.5, 6.5, 8.5, and 9.5. The relaxation is clearly exponential in all 
cases. Fig. 5 shows T c as a function of a^ 1 for simulated aggregates with 1642 
cells and surface tensions 4.5, 5.5, 6.5, 7.5, 8.5, and 9.5. A statistical analysis 
of the fit in the figure supports an intercept of the line with the vertical axis at 
the origin (p- value of linear coefficient = 0.8), consistent with Eq. ^ According 
to Eq. Q] the slope of the fit is the Rx viscosity, however simulations predict 
an exponent higher than 1 from a log-log plot of T c vs. aggregate radius. As 
we can see in Fig. 6, we determined that the slope for simulation is about 2.5 
and for chicken data is about 1.5 (see also fit in Fig. 2). Simulation exponent 
seems robust to variations in temperature and surface tension. In both cases 
the results suggest a deviation from hydrodynamics. 

3 Conclusions 

We have presented new experimental results on the relaxation of biological cell 
aggregates (chick embryonic neural retina aggregates). Our motivation was to 
measure the equivalent viscosity of this tissue. Such measurements are im- 
portant to characterize cell adhesiveness and to provide model parameters for 
theoretical or numerical models of collective motion ^7]. The relaxation of 3D 
aggregates is exponential, however, contrary to 2D aggregates of Hydra cells 
8 , relaxation time seems to increase faster than linearly with aggregate size. 
We are currently working on the verification of this experimental dependence 





6 



in case of 3D aggregates. To this end, we need to improve the experimental 
environment (CO2 and O2 control) in order to increase the observation time to a 
few days and be able to study larger aggregates. It should be also noted that 2D 
results suffer large errors. However if this dependence is confirmed, the reason 
why 2D Hydra aggregates behave as viscous fluid but not the 3D chick aggre- 
gates is not completely clear. The answer may arise from cell trajectory studies. 
In case of 2D hydra aggregates, all cells participated to rounding, motion was 
very collective |S] • We have not performed a similar study yet but it seems that 
in larger 3D aggregates, intense cell movements occur mainly in the aggregate 
periphery in a band of a few cell diameter while in the core of the aggregate cells 
are less active. Assuming nevertheless a fluid behavior, we obtain the following 
estimate for the tissue viscosity of 9 day old neural retina, 77 ~ (4 ± 2) x I0 6 
Pa.s. 

A simple Potts model simulation of the phenomena captures the exponential 
decay and indicates a superlinear dependence of relaxation time with aggregate 
size as observed in experiments with chicken cells. The simulation exponent, 
however does not agree with experiments but a possible explanation for the 
difference may come from cell trajectory studies. 3D simulations might also 
provide useful quantitative comparisons with 3D aggregates. 
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Figure 1: Shape relaxation in 2D Hydra aggregates (data from (A) 2D 

experimental set-up used to study aggregate shape relaxation. (B) A pure endo- 
dermal aggregate at Oh and (C) at 6h. Bar, 100 [im. (D) Deformation parameter 
d (see text) as a function of time (bullets: pure endodermal aggregate; triangles: 
mixed aggregate; circles: pure ectodermal aggregate). (E) Relaxation time T c 
as a function of aggregate radius (same legend as D). 



9 




Figure 2: Rounding of 3D chicken embryonic neural retina aggregates. (A) 
Configuration used to study shape relaxation: the 3D aggregates are floating 
near to the bottom of glass slide in culture medium at 37°C. (B) A small 9-day 
old neural retina aggregate (R=63 fim) at Oh and (C) at 24h. Bar, 50 nm. 
(D) Deformation parameter d as a function of time for three different aggregate 
sizes. (E) Relaxation time T c as a function of aggregate radius. Solid line is a 
power law fit (~ i? 152 , correlation coefficient, 0.939), dotted line is a linear fit 
with with Eq. [21 (correlation coefficient, 0.880). 
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Figure 4: Log-linear plot of deformation as a function of time for simulated 
aggregates with 1642 cells and surface tensions 4.5 (squares), 6.5 (crosses), 8.5 
(dots), and 9.5 (circles). 
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Figure 5: Plot of T c as a function of <r 1 for simulated aggregates with 1642 cells 
and surface tensions 4.5, 5.5, 6.5, 7.5, 8.5, and 9.5. The correlation coefficient 
of the fit is 0.85. 
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Figure 6: Log- log plot of T c vs. R for simulation and chicken data (triangles). 
R is lattice sites for simulation and fim for experiments. Simulation surface 
tensions are: 4.5 (circles), 5.5 (squares), 6.5 (crosses). An identity line is shown 
for slope comparison. 
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